NASA TN 0-7119 


NASA TECHNICAL NOTE 


^' 13-^33 



CASE FILE 

COPY 


ELASTOSTATIC STRESS ANALYSIS 
OF ORTHOTROPIC RECTANGULAR 
CENTER-CRACKED PLATES 

by George S. Gyekenyesi and Alexander Mendelson 

Lewis Research Center 
Cleveland, Ohio 44135 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • NOVEMBER 1972 



1. Report No. 2. Government Accession No. 

SA TN D-7119 


4. Title and Subtitle 

ELASTOSTATIC STRESS ANALYSIS OF ORTHOTROPIC 
RECTANGULAR CENTER- CRACKED PLATES . 


7. Author(s) 

George S. Gyekenyesi and Alexander Mendelson 


9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


3. Recipient's Catalog No. 


5. Report Date 
November 1972 


6. Performing Organization Code 


8. Performing Organization Report No. 

E-6889 


10. Work Unit No. 
501-22 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Note 


14. Sponsoring Agency Code 



16. Abstract 

A mapping- collocation method was developed for the elastostatic stress analysis of finite, 
anisotropic plates with centrally located traction- free cracks. The method essentially consists 
of mapping the crack into the unit circle and satisfying the crack boundary conditions exactly 
with the help of Muskhelishvili’s function extension concept. The conditions on the outer bound- 
ary are satisfied approximately by applying the method of least- squares boundary collocation. 

A parametric study of finite-plate stress intensity factors, employing this mapping- collocation 
method, is presented. It shows the effects of varying material properties, orientation angle, 
and crack- length-to-plate- width and plate- height- to- plate- width ratios for rectangular ortho- 
tropic plates under constant tensile and shear loads. 


17. Key Words (Suggested by Author(s)) 

Fracture mechanics 

Elasticity 

Anisotropy 


Cracks 

Plates 


18. Distribution Statement 

Unclassified - unlimited 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

32 


22. Price* 

$3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 





















ELASTOSTATIC STRESS ANALYSIS OF ORTHOTROPIC 
RECTANGULAR CENTER-CRACKED PLATES 
by George S. Gyekenyesi and Alexander Mendelson 
Lewis Research Center 

SUMMARY 

A mapping- collocation method was developed for the elastostatic stress analysis of 

v* , 

“finite, anisotropic plates with centrally located traction-free cracks. The method es- 
sentially consists of mapping the crack into the unit circle and satisfying the crack 
boundary conditions exactly with the help of Muskhelishvili's function extension concept. 
The conditions on the outer boundary are satisfied approximately by applying the method 
of least-squares boundary collocation. 

A parametric study of finite-plate stress intensity factors, employing this mapping- 
collocation method, is presented. It shows the effects of varying material properties, 
orientation angle, and crack- length- to- plate- width and plate- height-to-plate- width ratios 
for rectangular orthotropic plates under constant tensile and shear loads. 


INTRODUCTION 

Anisotropic materials are frequently applied in present-day structural design, as 
exemplified in composite materials. Therefore, the problems pertaining to the stress 
distribution in finite and quasi- infinite plates of such materials containing cracks are of 
immediate importance. The infinite- plate problems present no difficulty (refs. 1 to 8); 
however, the problems of finite anisotropic plates with cracks remained fairly neglected 
until quite recently (refs. 9 to 11). 

The mapping- collocation method as applied in this report was developed in refer- 
ence 10 and parallels Bowie’s mapping- collocation technique (ref. 12) for a finite iso- 
tropic plate with a central traction-free crack. 

This report analyzes the effects of varying material properties, orientation angle, 
and crack-length-to-plate-width and plate- height-to-plate- width ratios on both the 
opening- and sliding-mode stress intensity factors in a rectangular orthotropic plate 


with a centrally located traction- free crack under constant tensile and shear loads. 

In addition to the parametric investigation, a brief summary of the development of 
the mapping- collocation method is presented; however, for a detailed exposition the 
reader is referred to reference 10. 


SYMBOLS 


The following list contains the more commonly used symbols and their representa- 
tive meaning unless otherwise. defined in the text. In general, all symbols are defined 
when first introduced. 


A m 

iL 

n complex constant in Laurent series associated with positive ex- 
ponent terms 

a 

half- crack length 

B ln 

n th complex constant in Laurent series associated with negative ex- 
ponent terms 

C ij 

ij th element of elastic compliance matrix 

»ln 

n modified complex constant associated with positive exponent 
terms in Laurent series 

*2n 

n th modified complex constant associated with negative exponent 
terms in Laurent series 

*c 

complex constant depending on material properties 


modified complex constants depending on sum and difference of 
Ajj and respectively 

c 

half- width of rectangular plate 

E 11’ E 22 

Young's moduli of elasticity with respect to coordinate directions 

F ln’ F 2n’ G ln’ G 2n 

complex quantities utilized in construction of coefficient matrix in 
least- squares boundary collocation method 

G 12 

shear modulus of an orthotropic material 

Hi, Ha- 

nondimensional opening- and sliding-mode stress intensity factors , 
respectively ' 

il. 

half- height of a rectangular plate 

K P K 2 

opening- and sliding-mode stress intensity factors, respectively 


2 


M 


NN, NP 


x,y 

z 

z l’ z 2 

6 

^k 

^12 

M 1? M 2 

a , a , a 
xx’ yy’ xy 

<P 

<P i( z i)> <^ 2 (z 2 ) 


total number of collocation points on outer boundary of region with 
central crack 

truncation numbers on infinite sums associated with negative and posi- 
tive exponent terms in Laurent series, respectively 

constant shear stress applied to boundary of square plate 

constant tensile stress applied to boundary of rectangular plate 

rectangular coordinates 

conventional complex variable 

complex variables modified by material properties 

material orientation angle 

complex variables in parametric planes 

Poisson's ratio in an orthotropic material 

Lekhnitskii 's material parameters 

components of stress tensor, rectangular coordinates 

Airy's stress function 

complex stress potential functions 


BRIEF DESCRIPTION OF MAPPING-COLLOCATION METHOD 

Consider a doubly connected orthotropic region L defined by a crack T c and the 
region's outer boundary r, as shown in figure 1(a). The crack and its exterior can be 
mapped, as shown, into the unit circle y and its exterior by the mapping function 

J~ 2 2 

z + T z - a 

( 1 ) 

a 


where z is the usual complex variable, x + iy. 

We introduce two additional complex variables, z^ and Zg, defined by 


z k = 


1 - iM k 1 + iM k _ 

£z+ £z k = .1, 2 


where and 3X6 the roots of the characteristic equation 


( 2 ) 


3 


( 3 ) 


'22 


~ ^ 26 ^ + (^66 + ^ 12 ^ ~ ^ 16 ^ + 


'66 


' 12 ' 


'11 


M 4 = 0 


and the C- are the material compliances appearing in the generalized Hooke's law for 
an orthotropic material. 

The regions with the boundaries T ck and r k obtained by the transformation 
(2) can also be mapped into their respective ^-planes by the mapping functions 



k = 1,2 


( 4 ) 


as illustrated in figures 1(b) and (c). 




Figure 1. - Mapping of doubly connected regions L, Lj, and L 2 
into their corresponding parametric planes. 
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The Airy stress function <p (x, y) can be written in terms of two analytic functions 
of the complex variables z^ and z 2 - 

<p(x,y) = /?e[F 1 (z 1 ) + F 2 (z 2 )] (5) 

With the following definitions : 


<Pj( Zj) = 


dF ] 

dz* 


and 


cp 2 (z 2 ) = 


dF 5 

dz 0 


the zero traction conditions on the crack can be satisfied by applying Muskhelishvili's 
function extension concept across the unit circle (ref. 13). This is accomplished by 
expressing one of the unknown complex stress potentials in terms of the other. If 

^ 2 ^ 2 ) is c h° sen as 

Mo " Mi /i \ Mo ~ Mi 

+J -rn«2> < 6 > 

M2 “ M 2 ' 2 / M2 - M 2 


the crack boundary conditions will be satisfied for any choice of j), provided it is 
analytic. 

The substitution of expression (6) into the boundary conditions for the outer boundary 
of the region results in two conditions given completely in terms of <P -^y), ant * 

■^iUA 2 )- 

The determination of cp depends on its form of representation and the satisfac- 
tion of the boundary conditions on the outer boundary. Specifically, it is assumed that 
can be represented in the form of a truncated Laurent series, such that 


<P 


NP NN 

l«l>= A 10 + E A ln ? ? + Z B ln ? i n 

n=l n=l 


(7) 


5 


y 



Figure 2. -Tension problem of 
rectangular plate with central 
crack. 



Figure 3. - Shear problem of square plate 
with central crack. 


Because of the symmetry conditions for the orthotropic plate problems in question 
(figs. 2 and 3), n becomes necessarily odd and the final forms of the boundary condi- 
tions for the outer boundary are obtained as 


0 ^ #F n " *^*F 11 • • • ^ eF lNP “ ^ F 1NP • ‘ • ^ eF 2NN"^ F 2NN 
y ReG n - . . . ^G 1Np - -^» G 1Np • • • ^ eG 2NN " ^* G 2NN 




where NN and NP are odd. The terms in equations (8) are defined as 


6 


% = 


A 11 + B 11 
a(ii 2 - 


*in = 


*2n=‘ 


A 11 ' B 11 
*0*2 -T^) 

A ln 

a n (M 2 - TI 2 ) 

|ln 

a"(M 2 - H 2 ) 


* c = (M 2 - M 1 )(fi 2 ' M 2 )(M 2 - Mj) 


F 11 - ^2 


- H 2 ) ^z 2 - a 2 + (4 2 - Mj) V^- a2 + (4 2 ' **1> 


G 11 = U 


F ln = ^2 - 


F 2n = 0*2 


( 9 ) 


l(42-M2 ) V Z l- a2+ M 2^2- ^V z i- a2+ ^2- ^ 

M 2 )^i + y z 2 - a 2 ^ + (m 2 - Mj) ^ 2 + V 2 ! - * 2 ^ + 0*i ‘ **2^2 " V z 2 ' 

7I 2 )^i- V^) +(iI 2' M l> ^ z 2 • V z 2 • +(fi l‘ M 2^ z 2 + V z 2 ' 

G ln = Mi(ji 2 ‘ 1I 2 ) ^ z 1 + V Z 1 • a2 j + ^2 (?I 2 ' ^l^ z 2 + V z 2 " a2 j + 1I 2 (4 1 ' **2> ^2 * V Z 2 ‘ ^ j 

G 2 n = Mi(M 2 - M 2 ) - V Z 1 ‘ a2 j + M 2 (?I 2 ‘ M l } ^2 ' V z 2 " ^ j + ^2^1 - M 2> ^2 + V®2 " a ^j 


Equations (8) must now be solved for the unknowns, ^ n and ^n' ^he applicability of 
the least-squares boundary collocation method (refs. 14 and 15) to the problems at hand 
is readily obvious. Equations (8) are to be approximately satisfied at M points on the 
outer boundary, while the number of unknowns is given by NN + NP + 1. Then the re- 
sulting matrix equation can be written as 


7 



g <g = jr (10) 

2MX(NN+NP+1) (NN+NP+l)xl ' 2MX1 


For the least-squares boundary collocation method, the problem becomes 

g T g<g* = g T X (11) 

where ‘if* is the desired solution vector. 

The solution of equation (11) was carried out on a digital computer by declaring all 
quantities and operations in double precision and utilizing the Gauss-elimination tech- 
nique with complete pivoting. It should also be mentioned that the elements of the coef- 
ficient matrix g were scaled by the devisor 



After the approximate solution vector <g* is obtained, the computation of the stress 
intensity factors and the stress components is a matter of substitution of the elements 
of ‘if* into the following expressions : 


K t = — /?e 

1 VI 


NP NN 

<^2 - "l>( a *b + £ ran *ln-Z) nan *'2n| 

n=3,5,7 n=3 


K ? = - — Re | 

2 VI 


M 


NP NN 

2^2 ‘ M 1 ) ( a<r D + X nan< ^ln ‘ X nan ^2n j 

n=3, 5, 7 n=3 


( 12 ) 


(13) 
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n=3, 5, 7 
NN 



n=3, 5, 7 


9 



- 2/?e 



n =3,5,7 


r 





NN 



n=3, 5, 7 



ae) 


For a detailed discussion and complete derivation of the above method, refer to 
reference 10. 

RESULTS AND DISCUSSION OF SOLUTIONS OF TENSION AND SHEAR PROBLEMS 

Specifications of Parameters 

The solutions of the two problems shown in figures 2 and 3 are based on the follow- 
ing specifications: 

(1) The nondimensional stress intensity factors are defined as H k = K^/-y/a T t 
for tension and H k = K k /\/a T g for shear (k = 1, 2). 

(2) With the exception of the investigation of the effects of material properties on the 
stress intensity factors, the ratios 


E 


22 


E 


11 


0.366 


and 


10 



= 0.257 


were used throughout the entire analysis. These ratios correspond to the properties of 
fiberglass . 

(3) In studying the effects of varying material properties on the stress intensity fac- 
tors in both the tension and shear problems, the following parameter ranges were con- 
sidered: 


^22 

0.05 < — < 1 


0.001 < 


h = l 6 = 45° 
c 3 c 


(4) In studying the effects of varying plate aspect ratio h/c and crack- length-to- 
plate- width ratio a/c on the stress intensity factors in the tension problem, the param- 
eter ranges were chosen as follows: 


0 <- < 1 


0.25 < - < 133 


5 = 45 l 


(5) The effects of variation of the orientation angle 6 were also analyzed for both 
the tension and shear problems, considering the following parameter ranges: 


E-6889 


11 


0 <- s - 
c 3 


0° < 6 < 90° 



c 


Examination of Mapping-Collocation Method 

The available results which can be used to check the accuracy of this method are 
those for centrally cracked infinite plates (refs. 1 and 3) and the results of Kobayashi, 
Isida, and Sawyer (refs. 16 to 18) for finite isotropic rectangular plates. In order to 
gain a certain degree of confidence in the mapping- collocation method, both the results 
for the infinite orthotropic and finite isotropic plates were verified. 

The first verification is obtained by making the plate dimensions large compared to 
the crack length. 

This resulted in = 1.0000 and Hg = 0.0000. The corresponding analytical re- 
sults for an infinite plate are Hj = 1 and Hg = 0. It was, however, also noted that 
the satisfaction of the stress boundary conditions markedly improved with the use of 
more and more terms in the Laurent series expansion of (p j(Cj). Thus, the results of 
the infinite plates were verified for both the tension and shear problems. 

The verification of Kobayashi 's, Isida's, and Sawyer's opening- mode stress inten- 
sity factors for various finite isotropic rectangular plates in tension was carried out by 
setting the ratios 



and 


E 

E 


n 

22 


/ E ll 

\ 2G 12 



= 1 


12 


The value of the orientation angle had no influence on the results. The cases considered 
are shown in table I, where the various stress intensity factors are tabulated and can 
readily be compared with each other. As can be seen from table I, agreement with the 
known results is excellent. 

In order to economize the mapping- collocation method with regard to the use of the 
number of terms in the Laurent series expansion of two different truncation 

numbers NN and NP were assigned as limits of the sums in the various expressions. 
The numbers NN and NP are odd for orthotropic materials and correspond to the 
number of the negative and positive exponent terms in the Laurent series expansion of 
the stress potential function <p ^(^). Upon investigating various NN/NP ratios for the 
finite plate problems, it was found that a NN/NP ratio of 5/3 resulted in fairly fast con- 
vergence with rather well- satisfied boundary conditions. This NN/NP ratio was re- 
tained throughout the whole investigation of the problems given in figures 2 and 3. 

One of the most important indications of the degree of accuracy of the solutions ob- 
tained by the mapping- collocation method comes from the examination of the boundary 
stress between collocation points. Since the least-squares method of collocation results 
in the satisfaction of the boundary conditions only in the approximate sense, the magni- 
tude and frequency of oscillations of the boundary stresses serve as an indicator to the 
’’exactness" of the solution of the stress boundary- value problem. Thus, in each case, 
the boundary stresses were also examined. For example, in the case of a square plate 
in tension with a/c = 2/3 and 6 = 45°, the largest error in the boundary stresses was 
found to be about 4 percent and it occurred in a sudden manner at the corners of the 
plate. 

In addition, the pattern of convergence of the stress intensity factors Hj and H 2 
was also investigated. The results are shown in table II and figures 4 and 5. It was ob- 



Figure 4. - Typical pattern of convergence of opening-mode stress intensity factor for 
tensile loading. Plate size, h/c = 1; crack length, a/c = 1/2; orientation angle, 
6-45°. 
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Figure 5. -Typical pattern of convergence of sliding-mode stress intensity factor for 
tensile loading. Plate size, h/c = 1; crack length, a/c ■ 1/2; orientation angle, 

6 = 450 . 


served that the insensitiveness of H 1 and H 2 to a change in the number of terms of 
the Laurent series expansion of the stress potential is an excellent indicator to 

how well the boundary conditions are satisfied by the least- squares approximation. 

Upon considering various materials, plate sizes (aspect ratios), a/c ratios, and 
orientation angles under the loading conditions shown in figures 2 and 3, two conclusions 
became obvious regarding the convergence of the stress intensity factors: 

(1) Convergence of the stress intensity factors is definitely affected by the material 
parameters. It was found, for example, that for a square plate in tension with a/c = 
2/3 and 6 = 45°, the use of the dimensionless material constant 



with any ^ 22^ E 11 ratio resulted in oscillatory values of Hj and Hg. In this case, 
the number of unknowns was taken as 83, which gave excellent results for higher values 
of the dimensionless material constant. 
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(2) Convergence of the stress intensity factors is also affected by the a/c ratios. 
For a square plate in tension with 


E 

E 


n 

22 


/ E ll 

\ 2G 12 



0.257 


and with E 2 2 /Ejj = 0.366, a/c = 0. 9, and 6 = 45° with the use of 83 unknowns, it was 
found that the boundary conditions were rather badly satisfied. The method probably 
resulted in unconverged stress intensity factors. For smaller a/c ratios, well- 
converged values of and Hg were obtained. 


SOLUTION OF TENSION PROBLEMS OF ORTHOTROPIC 
RECTANGULAR PLATE WITH CENTRAL CRACK 

The consideration of the tension problem, that is, when the constant boundary 
stress is applied in the y-direction (fig. 2), resulted in various parametric studies in- 
volving material properties, a/c ratios, h/c ratios, and orientation angles. 


Effects of Material Properties 


An orthotropic material can be characterized by two dimensionless ratios con- 
structed from the engineering material constants E^p E 22 ? and (ref. 10). 

These dimensionless ratios are 


E 


E 


22 

11 


and 



The effects of these two dimensionless ratios on the stress intensity factors are shown 
in table HI and figures 6 and 7. 
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Figure 6. - Effects of material properties on opening-mode 
stress intensity factor for tensile loading, Plate size, 
h/c- 1; crack length, a/c“ 2/3; orientation angle, 
6=45°. 



Figure 7. - Effects of material properties on sliding-mode 
stress intensity factor for tensile loading. Plate size, 
h/c = 1; crack length, a/c» 2/3; orientation angle, 

6 * 45°. 


Figure 6 demonstrates the effects of the dimensionless material ratios on the 
opening-mode stress intensity factor, while figure 7 shows the sliding-mode stress in- 
tensity factor as a function of these two ratios. Figure 6 shows that, as the E 2 2 / E 1 1 
ratio decreases, the opening-mode dimensionless stress intensity factor Hj increases 
in value and that, for a fixed value of this ratio, the stress intensity factor exhibits a 
minimum value. 

Figure 7 shows the interesting result that there exists a sliding-mode stress inten- 
sity factor for the case of a pure tensile load, contrary to the results for an infinite or 
isotropic plate. This results strictly from the orthotropy of the material and the finite- 
ness of the plate dimensions. It may be noted that the maximum magnitude of the 
sliding-mode stress intensity factor shown in figure 7 is about 14 percent of the opening- 
mode stress intensity factor for the same E 22^11 ratio ' 


Effects of Plate Dimensions 


In order to investigate the effects of various plate widths and plate lengths on the 
stress intensity factors and Hg, the dimensionless material ratios were set at 


16 






= 0.257 




with an orientation angle 6 = 45°. These ratios correspond approximately to fiber- 
glass properties. The results are shown in table IV and figures 8 and 9. 


Effects of Orientation Angle and Crack Length 

For the investigation of the effects of the orientation angle and the crack length on 
the stress intensity factors, the h/c ratio was taken as 1 (square plate) and the dimen- 
sionless material ratios were held constant at the values specified previously. 

The variations of the orientation angle 6 and the a/c ratio resulted in table V. 
The graphs constructed from table V are presented in figures 10 and 11. Figure 10 
shows the effects of various orientation angles 5 on the opening-mode stress intensity 
factor Hj for constant a/c ratios. Figure 11 shows the sliding-mode stress intensity 


Crack length, 
a/c 



Orientation angle, 6, deg 


Figure 10. - Effects of orientation angle and plate width on opening- 
mode stress intensity factor for tensile loading. Plate size, h/c= 1. 
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Crack length, 



Figure 11. - Effects of orientation angle and plate width on 
sliding-mode stress intensity factor for tensile loading. 
Plate size, h/c = 1. 


factor H 2 as affected by changes in 6 and a/c. 

As an illustration of the typical stress distribution on the line y = 0 for the range 
1 < x/a ^ 3/2, the values of the stress components were tabulated in table VI for a 
plate with h/c = 1, a/c = 2/3, 6 = 45° and plotted in figures 12 and 13. 

The presence of a xy results from the general orthotropy of the material. It should 
also be observed that a xy becomes effectively zero on the boundary, which it should be 
for proper satisfaction of the boundary conditions. 

For each and H 2 value obtained, the complete boundary- value problem of an 
orthotropic rectangular plate containing a central crack had to be solved. In each case, 
in addition to the computation of the stress components on the boundary, the stress com- 
ponents on the y = 0 line were also recorded. It gives one confidence in the method to 
know that for isotropic and specially orthotropic materials, the shear stress component 
°xy was e ff ec tively zero on the y = 0 line. Also for problems with E 22 /E 11 = 1 , the 
shear stress component o_ on the y = 0 line was practically zero for any value of 

xy 


E 


E 


n 

22 


/ E ll 

V 2G 12 
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Figure 12. - Normal stress distribution on y ■= 0, x> a 
line for tensile loading. Plate size, h/c = 1; crack 
length, a/c = 2/3; orientation angle, 6= 45°; constant 
tensile stress, T t = 1. 



Figure 13. - Shear stress distribution on y=0, x>a line 
for tensile loading. Plate size, h/c = 1; crack length, 
a/c = 2/3- orientation angle, 6= 45°; constant tensile 
stress, Tj » 1. 


SOLUTION OF SHEAR PROBLEM OF ORTHOTROPIC 
SQUARE PLATE WITH CENTRAL CRACK 

In addition to the tension problem discussed in the preceding section, the problem of 
an orthotropic square plate containing a central crack and loaded by unit shear stress 
(fig. 3) was also considered. As in the case of the tension problem, the shear problem 
was also solved for various dimensionless material ratios, a/c ratios, and orientation 
angles. The h/c ratio was held constant (h/c = 1) throughout the entire analysis of the 
shear problem. 

The results of the parametric study involving variations of the dimensionless mater- 
ial ratios while we held a/c - 2/3 and 6 = 45° are summarized in table VII and fig- 
ures 14 and 15. 
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.01 .1 1 10 100 


E 11 ,E 22 



Figure 14. - Effects of material properties on opening-mode 
stress intensity factor for shear loading. Plate size, 
h/c = 1; crack length, a/c° 2/3; orientation angle, 8 = 
45°. 



Figure 15. - Effects of material properties on sliding-mode 
stress intensity factor for shear loading. Plate size, 
h/c= 1; crack length, a/c° 2/3; orientation angle, 6' 
45°. 


The effects of variation of the orientation angle 6 and the a/c ratio are shown in 
table VTTT and figures 16 and 17. There is an important observation which concerns the 
effects of the variation of the orientation angle in the shear problem. As far as 5- 
dependence is concerned, the orthotropic plate behaves differently in tension and in 
shear. In the tension problem it was found that the maximum Hj values always oc- 
curred at 6 = 90° (specially orthotropic material), while in the shear problem it is 
obvious that the maximum H 2 is 6-dependent. However, in both the tension and shear 
problems the minimum values of the significant stress intensity factors occurred at 
6 = 0°, which designates the other type of special orthotropy. 

As an illustration of the typical stress distribution on the line y = 0, 1 < x/a 2 = 3/2, 
the values of the stress components were recorded in table EX and depicted in figures 18 
and 19 for a plate with a/c = 2/3 and 6 =45°. Both and cr yy were found to be 
zero on the line y = 0, 1 < x/a < 3/2 with ^ 22^11 ~ 1 and for any value of 
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Normal stress 




The special cases of isotropy and special orthotropy resulted also in zero values for the 
normal stresses and a ^ on the y = 0, 1 < x/a s 3/2 line. 

CONCLUDING REMARKS 

The mapping- collocation method as developed in reference 10 was applied to a large 
number of rectangular orthotropic plate problems in order to study the effects of vary- 
ing material properties, orientation angles, and crack- length-to- plate- width a/c and 
plate-height- to- plate- width h/c ratios on the stress intensity factors for both tensile 
and shear loadings. 

A natural extension of the mapping- collocation method as applied to finite ortho- 
tropic regions with centrally located traction- free cracks would be to consider various 
shapes for the inner boundary. For example, triangular, rectangular, and elliptical 
boundaries could be specified by a single general Schwarz- Christoffel transformation 
and mapped into the unit circle. This possibility would require further research as to 
the accuracy of the method when it is applied to various types of doubly connected re- 
gions. 

Another direction of research presents itself in the consideration of a more detailed 
parametric study of the effects of dimensionless material constants, orientation angles, 
a/c ratios, and h/c ratios on the stress intensity factors. It is a possibility that, for 
certain parameter ranges, "practical forms" of expressions could be obtained for the 
approximation of the stress intensity factors. These forms would result from curve 
fitting the various parametric data. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 21, 1972, 

501-22. 
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TABLE I. - COMPARISON OF OPENING-MODE STRESS INTENSITY 


FACTORS FOR RECTANGULAR ISOTROPIC PLATES 
WITH CENTRAL CRACKS 


[Plate size, h/ c = 2 7) 


Crack length, 
a/c 

Opening-mode nondimensional stress intensity factor, Hj 

Kobayashi 
(ref. 16) 

Isida 
(ref. 17) 

Sawyer 
(ref. 18) 

Present 

study 

1/12 

1.0048 

1.0040 

1.0094 

1.0042 

1/6 

1.0191 

1.0154 

1.0253 

1.0171 

1/4 

1. 0448 

1.0392 

1.0509 

1.0395 

1/3 

1.0830 

1.0750 

1.0737 

1.0733 

1/2 

1.2215 

1. 1725 

1.1443 

1.1877 

2/3 

1.4665 

1.4142 

1.3804 

1.4162 


TABLE II. - TYPICAL PATTERN OF CONVERGENCE OF STRESS INTENSITY 
FACTORS FOR TENSILE LOADING 


[Plate size, h/c = 1; crack length, a/c = 1/2; orientation angle, 6 = 45°!] 


Number of 
equations 

Number of unknowns 

15 

25 

41 

83 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

16 

24 

32 

40 

48 

64 

96 

128 

176 

1.3548 

1.3466 

1.3455 

-0.0377 

-.0394 

-.0397 















N 




1.3531 

-0.0461 







1.3567 

-0.0475 



1.3451 

- . 0400 

1.3520 

- . 0460 



1.3562 

1.3562 

-.0472 
- . 0472 



1.3451 

-.0402 

1.3520 

-.0461 

1.3563 



-0. 0472 










TABLE m. - EFFECTS OF MATERIAL PROPERTIES ON STRESS INTENSITY 
FACTORS FOR TENSILE LOADING 


[Plate size, h/c = 1; crack length, a/c = 2/3; orientation angle, 6 = 45°|] 


E ll/ E 22 

E 22' /E ll 

/ e a \ 2 

0.05 

0.2 

0 

.4 

0.6 

0.8 


1 

.0 

















V 2G 12 7 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 

2 

0. 001 

2.475 

-0. 651 

2.358 

-0.345 

2.339 

-0. 189 

2.320 

-0.102 

2.318 

-0.006 

2.392 

-0. 

032 

.050 

2.006 

-.387 

1.805 

-.184 

1.741 

-.099 

1.735 

-.057 

1.731 

-.024 

1.730 

0 


.250 

1.963 

-.331 

1.726 

-.153 

1.663 

-.084 

1. 644 

-.045 

1.636 

-.020 

1.634 



.500 

1.964 

-.318 

1.712 

-.144 

1.646 

-.078 

1.624 

-.042 

1.617 

-.018 

1.614 



. 750 

1.967 

-.311 

1.709 

-.140 

1.639 

-.076 

1.618 

-.041 

1.610 

-.017 

1.607 



1.000 

1.972 

-.307 

1.707 

-.139 

1.636 

-.075 

1.614 

-.041 

1.605 



1.604 



2. 000 

1.983 

-.300 

1.706 

-.134 

1.634 

-.072 

1.610 

-.040 

1.601 



1.600 



4. 000 

1.994 

-.296 

1.710 

-.131 

1.635 

-.071 

1.610 

-.040 

1.601 



1. 602 



10. 000 

2.008 

-.291 

1.714 

-. 129 

1.638 

-.069 

1.614 

-.038 

1.604 



1.603 



100. 000 

2.034 

-.287 

1.727 

-.127 

1.648 

-.069 

1.622 

-.038 

1.612 



1.611 



CO 

2.048 

-.284 

1.735 

-.127 

1.655 

-.069 

1.629 

-.038 

1.619 



1.618 




TABLE IV. - EFFECTS OF CRACK LENGTH AND PLATE SIZE ON STRESS INTENSITY 

FACTORS FOR TENSILE LOADING 


[Orientation angle, 5 = 45 °J 


Crack 





Plate size 

h/c 







length, 

a/c 

1/4 

1/2 

1 

2 

4 

133 


H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

H 1 

H 2 

0 

1.000 

0.000 

1.000 

0.000 

1.000 

0.000 

1.000 

0. 000 

1.000 

0.000 

1.000 

0 

1/6 

1.417 

-.055 

1. 117 

-.016 

1.040 

-.006 

1.023 

-.006 

1.023 

-.006 

1. 006 



1/3 

2.331 

-. 199 

1.450 

-.062 

1.157 

-.021 

1.093 

-.021 

1. 089 

-.021 

1. 017 



1/2 

3.465 

-.400 

2.020 

-.151 

1.356 

-.047 

1.232 

-.045 

1.214 

-.045 

1.038 



2/3 

4. 946 

-.680 

2.999 

-.288 

1.669 

-.091 

1. 497 

-.082 

1.460 

-.076 

1.069 
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TABLE EX. - STRESS DISTRIBUTION ON y = 0, x > a 
LINE FOR SHEAR LOADING 


[Plate size, h/c = 1; crack length, a/c = 2/3; orientation 
angle, 6 = 45°; constant shear stress, T g = ll] 


x/a 

Components of stress tensor (rectangular coordinates) 


0 

1 

O 

0 


XX 

yy 

xy 

1.016 

-3.01 

-0. 60 

8.77 

1.046 

-1.55 

-.36 

5.14 

1.076 

-1.09 

-.28 

4.02 

1.106 

-.085 

-.23 

3.43 

1. 136 

-.68 

-.19 

3.04 

1.166 

-.56 

-.15 

2.76 

1.196 

-.46 

»— ( 
rH 

1 

2.53 

1.228 

-.37 

-.07 

2.34 

1.258 

-.29 

r 

O 

CO 

2. 17 

1.288 

-.23 

.03 

2.01 

1.318 

-.18 

.09 

1.87 

1.348 

-.12 

. 17 

1.72 

1.374 

-.08 

.25 

1.58 

1.410 

-.04 

.35 

1.44 

1.440 

-.02 

.46 

1.30 

1.470 

0 

.58 

1. 15 

1.500 

0 

.72 

1.00 
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